home *** CD-ROM | disk | FTP | other *** search
/ Libris Britannia 4 / science library(b).zip / science library(b) / MATH / FSULTRA1.ZIP / ULTRA_TP.ASM < prev    next >
Assembly Source File  |  1992-06-18  |  3KB  |  143 lines
  1. comment ! 
  2. FSU - ULTRA    The greatest random number generator that ever was
  3.         or ever will be.  Way beyond Super-Duper.
  4.         (Just kidding, but we think its a good one.)
  5.  
  6. Authors:    Arif Zaman (arif@stat.fsu.edu) and
  7.         George Marsaglia (geo@stat.fsu.edu).
  8.  
  9. Date:        27 May 1992
  10.  
  11. Version:    1.05
  12.  
  13. Copyright:    To obtain permission to incorporate this program into
  14.         any commercial product, please contact the authors at
  15.         the e-mail address given above or at
  16.  
  17.         Department of Statistics and
  18.         Supercomputer Computations Research Institute
  19.         Florida State University
  20.         Tallahassee, FL 32306.
  21.  
  22. See Also:    README        for a brief description
  23.         ULTRA.DOC    for a detailed description
  24.  
  25. -----------------------------------------------------------------------
  26. ! 
  27. ;
  28. ; File:  ULTRA_TP.ASM (Turbo Pascal calling conventions)
  29. ;
  30. DOSSEG
  31. .MODEL TPASCAL
  32.  
  33. ;===== B: GLOBAL FUNCTIONS, SUBROUTINES and VARIABLES =================
  34. ;
  35.     public  i31bit, i15bit, i7bit,  i1bit,  uni,    duni
  36.     public  i32bit, i16bit, i8bit,  rinit,  vni,    dvni
  37.  
  38. ;===== D. MACRO DEFINITIONS ===========================================
  39. ;
  40. ; RinitProcStart should take two 32-bit arguments conx and shrx
  41. ;   and place them in the data segment in congx and shrgx.
  42. ;   es and ds should both point to the data segment.
  43. ;   conx must be odd (so we or with 1 just to make sure).
  44. ; FillProc
  45. ;   DS is already the data segment. ES should also be made the same.
  46.  
  47. EnterProcedure   macro
  48. endm
  49.  
  50. ExitProcedure    macro
  51.     ret
  52. endm
  53.  
  54. Enter2arg   macro
  55.   arg conx:dword, shrx:dword
  56.   EnterProcedure
  57.     mov ax,ds
  58.     mov es,ax
  59.     mov ax,word ptr conx
  60.     shl ax,1
  61.     mov conglo,ax
  62.     mov ax,word ptr conx+2
  63.     rcl ax,1
  64.     mov conghi,ax
  65.     or  byte ptr conglo,1
  66.     mov ax,word ptr shrx
  67.     mov shrglo,ax
  68.     mov ax,word ptr shrx+2
  69.     mov shrghi,ax
  70. endm
  71.  
  72. Exit2arg     macro
  73.   ExitProcedure
  74. endm
  75.  
  76. EnterFill    macro
  77.     mov    ax,ds
  78.     mov    es,ax
  79. endm
  80.  
  81. ExitFill    macro
  82.     ret
  83. endm
  84.  
  85. DwordFn macro
  86. endm
  87.  
  88. WordFn  macro
  89. endm
  90.  
  91. ByteFn  macro
  92. endm
  93.  
  94. RealFn  macro
  95. endm
  96.  
  97. DoubleFn macro
  98. endm
  99.  
  100. ;===== DATA ========================================================
  101.  
  102.   N equ 37              ; The number of 32 bit words in the table
  103.   N2 equ 24             ;  x(i) = x(i-N) - x(i-N2)
  104.  
  105. counterpoint    STRUC
  106.   c  dw 0               ; counter
  107.   p  dw ?               ; offset for pointer
  108.      dw ?               ; segment for pointer
  109.   x  dd N dup (?)    ; data
  110. ENDS
  111.  
  112.   extrn swbstate            ; In TP, the shared DATA is
  113.   cps = size counterpoint
  114.   cp1 = cps - 4*N + 32
  115.   swb32 = counterpoint ptr swbstate    ; in the TPU, so we define a
  116.   swb16 = counterpoint ptr swb32+cps    ; large array there, and refer
  117.   swb8  = counterpoint ptr swb16+cps    ; to it here. This ensures that
  118.   swb1  = counterpoint ptr swb8 +cps    ; the swbstate is contiguous.
  119.   swbseed = dword ptr swb1     + cp1
  120.   flags   = byte  ptr swbseed  + 4*N
  121.   congx   = dword ptr flags    + 1
  122.  
  123.   extrn tmpq  : qword
  124.   extrn neg31 : word
  125.   extrn neg63 : word
  126.  
  127.   conglo  equ word ptr congx
  128.   conghi  equ word ptr congx+2
  129.   tmpw1   equ word ptr tmpq
  130.   tmpw2   equ word ptr tmpq+2
  131.   tmpw3   equ word ptr tmpq+4
  132.   tmpw4   equ word ptr tmpq+6
  133.   tmpdlo  equ dword ptr tmpq
  134.   tmpdhi  equ dword ptr tmpq+4
  135.   shrglo  equ tmpw1
  136.   shrghi  equ tmpw2
  137.   shrgx   equ tmpdlo
  138.  
  139. .CODE
  140.   INCLUDE ULTRACOD.INC
  141. END
  142.  
  143.